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Abstract 

The leading /?Q _1 a" terms in the Wilson coefficient and anomalous dimension of 
the chromomagnetic operator in the heavy-quark effective Lagrangian are summed 
to all orders of perturbation theory. The perturbation series for the anomalous di- 
mension is well behaved, while that for the Wilson coefficient exhibits a divergent 
behaviour already in low orders, caused by a nearby infrared renormalon singu- 
larity. The resulting ambiguity is commensurate with terms of order 1 / m 2 in the 
effective Lagrangian, whose corresponding ultraviolet renormalons are identified. 
An excellent approximation for the scheme-invariant Wilson coefficient at next- 
to-next-to-leading order in renormalization-group improved perturbation theory is 
obtained. 
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1 Introduction 



The properties of hadronic bound states containing heavy quarks are characterized by 
a large separation of energy scales. Effects associated with the heavy-quark mass m 
are perturbative and can be controlled once they have been separated from other, long- 
distance effects. This separation is most conveniently done using an effective low-energy 
theory. In the case of hadrons containing a single heavy quark, the relevant effective 
theory is the heavy-quark effective theory (HQET) (see |jj for a review). Its Lagrangian 
is 

h v iv-Dh v + J- h v (iD ± ) 2 h v + \ ' ^ h v g s a^ lu h v + 0(l/m 2 ) , (1) 



2m 4m 
where v is the velocity of the hadron containing the heavy quark, h v the velocity- 
dependent heavy-quark field, and D± the covariant derivative orthogonal to v. The 
operators arising at order 1/m in the effective Lagrangian correspond to the kinetic en- 
ergy of the heavy quark and its chromomagnetic interaction 0, [|. Hadronic matrix 
elements of these operators appear in many applications of HQET. The coefficient of 
the kinetic-energy operator is fixed by reparametrization invariance, an invariance un- 
der infinitesimal changes of the velocity Hence, the only non-trivial short-distance 
coefficient in the effective Lagrangian at next-to-leading order in 1/m is the coefficient 
C (m, /i) of the chromomagnetic operator. 

The dependence of the Wilson coefficient on the large scale m can be factorized using 
the renormalization group. In general, 



a s (m) 

C(m,ri=C(m, m)exp / da 
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is the /3-function, and 
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(2) 



(3) 
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the anomalous dimension of the chromomagnetic operator (given in the MS scheme). 
The two-loop coefficient 71 of the anomalous dimension was calculated in || [|, and the 
two-loop initial condition C(m,m) in 0. The corresponding one-loop expressions were 
obtained much earlier in ||. The factor K{jj) in (^|) compensates the scheme and scale 
dependence of the chromomagnetic operator in the effective Lagrangian (|I|). All short- 
distance physics relevant to physical observables is contained in the renormalization- 
scheme invariant function C(m) which, up to an m- independent normalization, is given 
by 



C(m) 
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In the MS scheme, the next-to-leading correction is given by H 
2 V S -!- 2(',-j - [ — C A 4 3( ',) - 
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f ^ + 3C F ) ^ + (7CU + 11C F ) ^ 



(6) 



In the above equations, Cp = MN — l/N), Ca = N, and Tp = 1/2 are the colour factors 
for an SU(iV) gauge group. 

As an application, consider the mass splitting between the ground-state pseudoscalar 
and vector mesons containing a heavy quark. Its leading contribution comes from the 
chromomagnetic operator in the effective Lagrangian, which is the first term breaking 
the heavy-quark spin symmetry. Including terms of order 1/m 2 , one obtains [[?]][] 
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p G (p) 



C(m, p) pl G {n) + C 2 (m, p) p\{p) - C LS (m, p) p\ s {p) 



(7) 



where p G {p) is the matrix element of the chromomagnetic operator between ground-state 
mesons. Similarly, the hadronic parameters pf (p) are defined in terms of the matrix el- 
ements of operators appearing at order 1/m 2 in the effective Lagrangian (7j: ppsil 1 ) 
parametrizes the spin-orbit interaction, while p^ G (p) and p\(p) are bilocal matrix ele- 
ments of kinetic-chromomagnetic and chromomagnetic-chromomagnetic insertions. (We 
omit the contributions from 4-quark operators ||, which are irrelevant to our discussion.) 
These parameters are independent of the heavy-quark mass, but they depend on the 
renormalization scale in such a way that the scale dependence of the Wilson coefficients 
is compensated. The coefficient of the spin-orbit term is related by reparametrization 
invariance to that of the chromomagnetic operator [ID -|]12| (see also the appendix of 
||) : Cps{ m ip) — 2C(m,p) — 1. Introducing then the renormalization-scheme invariant 
parameters 



we find 
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All short-distance effects in this relation are contained in the single coefficient C(m), 
which is the object of our study. 



1 A similar mass formula has been derived in [g]; however, the term proportional to p\ is missing 
there. 
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2 All-order results in the large-/3o limit 



The perturbation series for C(m, fi) can be arranged as 

oo L— 1 

C(m,fi) = 1 + E E aPWriPS «* l (a0 , (io) 

l=i n =o 

where L is the number of loops, and /3q the leading coefficient of the /5-function in (|3|). 
In this letter, we sum the terms of order Pq -1 ^ to all orders of perturbation theory. In 
other words, we consider the limit of large @ for fixed p a s and calculate the coefficient 
C(m, /i) to order l//3 , neglecting terms of order 1/Pq and higher. Strictly speaking, 
there is no sensible limit of QCD in which f3 may be considered a large parameter 
(except, maybe, rif — > — oo); however, retaining only the leading (3q terms often gives 
a good approximation to exact multi-loop results (see, e.g., Hl3|), in particular in cases 
when there is a nearby infrared renormalon |TJ]]. At the least, it will provide us with 



some information about the summability of the perturbation series. 

The coefficients a^li of the terms with the highest degree of /3q in (10) are determined 



by diagrams with L — 1 light-quark loops, which are rather straightforward to calculate. 
We work in dimensional regularization with d = 4 — 2e space-time dimensions and adopt 
the MS subtraction scheme. At first order in 1/Pq, multiplicative renormalization simply 
amounts to a subtraction of the l/e n poles, and coupling-constant renormalization is 
given by (ft 2 = /i 2 e 7 /47r) 

M _ -2e b _ PqoJjj) _ 1 ri1 x 

(4tt)2 A* 1 + b/£ , 4n 21n(A*/A H! j)- 1 J 

The perturbation series for C(m,/i) can then be written as 

1 00 Fie Le) ( b \ L 
C(m, fj.) = 1 + — L ~ ( minimal subtractions) + 0(1/ ffi) . (12) 



The function F(e, u) is regular at e = u = 0. Following the methods of |1| [0J (used also 



in jl3|), we now expand F(e,u) in powers of e and u, and [b/(e + b)] L in powers of b/e, to 
obtain a quadruple sum in (p!2|). Combinatoric identities relate the 1/e terms, and hence 
the anomalous dimension of the chromomagnetic operator, to the Taylor coefficients of 
F(£,0)|T§: 

2b 

7 = -F(-6,0)+O(l//3 2 ). (13) 

Po 

The finite terms, which determine the Wilson coefficient itself, receive contributions from 
the Taylor coefficients of both F(e, 0) and F(0, u) ||16|| : 
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In this expression, all dependence on the heavy-quark mass resides in the function 
F(0,u) ~ (fi/m) 2u , while a dependence on the renormalization scale also enters through 
the coupling b ~ a s (fi). By separating in (|I4T) terms depending on the two scales m and 
fi, we find that the next-to-leading logarithmic correction c(m) in (|5|) is given by 



1 

:(m) = -L [ due- u[1/b{m)+K] S(u) + 0(1/ p 2 ) , 



F(0,u)-F(0,0) 



u 



(15) 



fi=m 



The same result can also be derived by evaluating (Q) in the large-/?o limit. Here 6(m) = 
f3oa s (m)/4:7r, and the constant k is introduced to compensate the scheme dependence 
of l/b(m). The most natural choice is to have k = —5/3 in the MS scheme. Then the 
combination 

1 + k = 2 \n(m/A m ) - - = 2 ln(m/A v ) (16) 



b(m) 



e 5 / 6 A m . The 



defines the inverse coupling in the so-called V scheme |[L7|| , with Ay 
function S'(-u) is the Borel transform of the perturbation series for c(m) in that scheme. 

The coefficient C(m,fi) is obtained by matching scattering amplitudes of an on-shell 
heavy quark in an external field in QCD and HQET, including terms of order 1/m @, 
As mentioned above, our focus is on L-loop diagrams with L — 1 light-quark loops. All 
HQET diagrams vanish because they contain no mass scale. The relevant QCD diagrams 
are shown in Fig. |l|. They must be supplemented by the wave-function renormalization of 
the external quark fields [0 . We find that in the MS scheme the result for the function 
F(e, u) in (O) has the form 



F(e,u) = (£) 2 V r(1 T ,| M)r(1 ] U) D(s)^- 1 \C F N F (e,u) + C A N A (s,u)} , (17) 



where 



r(3-tt-e) 



1 ' r(4-2e) 3 + 



(18) 



is related to the contribution of a light-quark loop to the gluon self-energy, and 

N F (e, u) = 4w(l + u- 2eu) , 
2-u-e r 



N A (e,u) 



2(1 



(2-5e + 2e 2 ) + (3 - Qe + Ae 2 )u 



(19) 



These formulae reproduce the known L = 1 and L = 2 || results, where L = u/e. 

Having obtained an explicit result for the function F(e, u), we are ready to derive all- 
order results for the anomalous dimension and Wilson coefficient of the chromomagnetic 
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Figure 1: QCD diagrams for heavy-quark scattering in an external field. A 
gray circle with a number n inside represents a chain on n light-quark loop 
insertions on a gluon propagator. 



operator in the large-/?o limit. From (0), we obtain 
2C A b(l + 26)r(5 + 2b) 
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(20) 



which reproduces correctly the leading (in /3 ) term of the exact two-loop result (^). The 
radius of convergence of the perturbation series in ( p0[) is /?o | ck s | < 47r. 

Next, setting // = m in Q14D , we find that the perturbative expansion of the Wilson 
coefficient at the scale \x = m reads 



C(m, m) = 
V317 



I f ^{ 2 (C F + C. 4 ) 



47T 



25 r + [ 299 -u 7 Hr 
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(21) 



which again reproduces the leading term of the exact two- loop result ||. Unlike the 
case of the anomalous dimension, here retaining only the leading term in 1/(3q gives a 
reasonable approximation; for iV = 3 colours, the exact two-loop result is 



C(m, m) 



1 + H + (2.869^o - 9.761) 
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7T 



(22) 



showing that the sub-leading term is about one third of the leading one. 

Finally, we quote the result for the next-to-leading logarithmic correction c(m) in @ 



Expanding the integral in (15) in powers of the coupling constant, we find the expansion 
coefficients 

25 

Cl = 2C F + -C A + O(l/(3 ), 
6 
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(23) 



The result for c\ reproduces the leading term in (|6]). For A/" = 3 colours, the numerical 
values of the coefficients are c x « 15.167, c 2 w 45.147/3 , c 3 « 226.64/3^, c 4 w 1526.0/$, 
C5 ~ 13035/5q, etc. They exhibit a (same-sign) factorial divergence c n+ i ~ (2j3 ) n n\, 
which renders the perturbation series for c(m) divergent and not Borel summable. On 
the contrary, the coefficients in the expansion of the anomalous dimension in ( p0|) appear 
to be well behaved; in fact, this series is rapidly convergent. Although we do not expect 
that the higher-order coefficients obtained in the large-/?o limit accurately represent the 
exact values, we do trust in the qualitative features of the results. This allows us to derive 
an approximation for the next-to- next-to-leading order (NNLO) correction C2 in (|J). 
Expanding (Q) to second order in a s (m), and inserting the known expansion coefficients 
for the three-loop /3-function the two- loop anomalous dimension || |6j, and the two- 
loop matching coefficient C(m,m) ||, we obtain an exact expression for c 2 depending 
only on the unknown three-loop coefficient of the anomalous dimension, for which we 
write 72 = Ca(— 0q + ?7iA) + Vo), where 771 and r] are unknown. This parametrization 
takes into account that the anomalous dimension vanishes in the abelian limit Ca = 
0]. We then expand c 2 = d_ij3 + J2 n =o ^n(CU//?o) n and find 
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2- 18 ■ (24 » 

where the numerical values refer to AT = 3 colours. Because of the good convergence 
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of the perturbation series for the anomalous dimension indicated by our analysis of the 
large-/3 limit, it is conservative to assume that \rjiP + r] \ < /3q, which is equivalent to 
the statement that the true value of the three-loop anomalous dimension differs from its 
value in the large-/?o limit by less than 100%. Under this assumption, we obtain 

03(71/ = 4) = 210.08 ±6.25, c 2 {n f = 3) = 238.04 ± 6.75 . (25) 

The uncertainty due to the unknown terms in the three-loop anomalous dimension is 
negligible compared with the overall size of the coefficient. Using this result together 
with the exact one- loop coefficient C\ in (((]), we find for 71/ = 4 

C(m) = [a s (m)} 9/25 [l + 0.672a s (m) + (1.33 ± 0.04) a* (m) + . . .] . (26) 

In the context of the heavy-quark expansion, this is the most precisely known Wilson 
coefficient to date. 



3 Borel summation and renormalon ambiguities 

As a consequence of its divergent behaviour, the perturbation series for c(m) must be 
truncated close to its minimal term, and the perturbative result for the Wilson coef- 
ficient is intrinsically ambiguous. Although we have explored this feature only in the 
large-/3 limit, it is believed to be of a rather general nature (see, e.g., In regular- 

ization schemes without an explicit infrared cutoff (such as dimensional regularization 
with minimal subtraction), the perturbative calculation of the coefficients c n involves 
an integration over all gluon momenta, including long-distance contributions from soft 
gluons. High orders in the expansion probe the region of increasingly smaller gluon 
momenta, a regime where perturbation theory is bound to fail. Remarkably, the pertur- 
bation series knows about its deficiency and signals it through the divergent behaviour 
of the expansion coefficients. We will now investigate this phenomenon in more detail. 

A convenient tool to study the ambiguities encountered in the attempt to resum an 
asymptotic perturbation series is provided by the Borel image of that series, given by 
the function S(u) in the integral representation fll5|). Using fll?D , we find that 



T(u)T(l-2u) 
S{U) = T(3-u) 



4«(1 + u)C F + ^ '-C/ 



- e" 5 "/ 3 ^ . (27) 
u 



If it existed, the Laplace integral in fllBD would define the Borel sum of the perturbation 
series for c(m) in the large-/3o limit. However, the presence of singularities along the 
integration contour (i.e. for positive values of u) makes the integral ill-defined. In the 
large-/3 limit, the function S(u) has pole singularities at half-integer values of u called 
infrared renormalons |2(| |2l|. Any attempt to resum perturbation theory involves an 



arbitrary choice of how to deal with these singularities. A measure of the ambiguity in 
the value of the Borel sum, which is of the same size as the minimal term in the series, is 
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provided by the residue of the nearest singularity, in the present case located at u — 1/2. 
Using (|i~6l), we obtain 

A[c W ] = 1(2 Cf + I C ,)^, (28) 
with Ay = e 5 / 6 Aj^g. The situation is similar to the well-known case of the infrared renor- 



malon ambiguity in the definition of the pole mass of a heavy quark p2| , |23| . The pertur- 
bation series relating the pole mass to a short-distance mass is also affected by an infrared 
renormalon at u — 1/2, and the corresponding ambiguity is Am = 2CpAy / /3q. Taking 
into account that the chromomagnetic operator appears in the effective Lagrangian (|I|) 
multiplied by a power of the inverse pole mass, we find that 



A 



C(m) 



m 



7( ' \ Am f 7 1 i , 2 



8C F m 2 



+ 0(l/ft). (29) 



Hence, the infrared renormalon ambiguity in the product C(m)/m is purely non-abelian, 
and commensurate with the contributions of higher-dimensional operators in the 1/m 
expansion. Indeed, in all predictions for physical quantities, the infrared renormalon 
ambiguities must cancel against corresponding ambiguities in the long-distance matrix 
elements of some higher- dimensional operators. 

When the heavy-quark expansion is applied to calculate a physical quantity, e.g. the 
mass splitting in (|7|), the resulting expressions contain short- distance Wilson coefficients 
and long-distance hadronic matrix elements. As mentioned above, in regularization 
schemes without a hard momentum cutoff the Wilson coefficients also contain contribu- 
tions from large distances, where perturbation theory is ill-defined, and these contribu- 
tions produce infrared renormalon ambiguities. Likewise, the hadronic matrix elements 
(which are not calculable perturbatively) contain contributions from small distances, 
which lead to ultraviolet renormalon ambiguities. In other words, in such schemes the 
separation of short- and long-distance contributions into Wilson coefficients and matrix 
elements is intrinsically ambiguous. Only when all contributions are combined to form a 
physical quantity, an unambiguous result is obtained. In the context of the HQET, the 
cancelations between infrared and ultraviolet renormalon ambiguities have been traced 



in detail in []24j. Here we consider the particular case of the mass splitting. 

The ultraviolet contributions to the hadronic parameters pf are independent of the 
nature of the external states and thus may be calculated using quark and gluon rather 
than hadron states. At one-loop order, the relevant diagrams are shown in Fig. [2| On 
dimensional grounds, they are linearly divergent in the ultraviolet region. In dimensional 
regularization such power divergences are not seen (by definition), but they reflect them- 
selves in a factorial divergence of the perturbation series when higher-order diagrams are 
taken into account. We shall again analyse these contributions in the large-/9o limit, by 
inserting an arbitrary number of light-quark loops in all possible ways into the diagrams 
of Fig. [| The relevant matrix elements are then expanded to linear order in the gluon 
momentum q and projected onto the structure of the chromomagnetic operator. The 
external quarks are kept off-shell in order to provide for an infrared regulator. The Borel 
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(al) (a2) (a3) (a4) 




(cl) (c2) (c3) 




(dl) 



Figure 2: Diagrams contributing to the ultraviolet renormalons contained in 
the parameters p\ G (a), p\ (b), and p\ s (c). Insertions of the chromomagnetic, 
kinetic and spin-orbit operators are indicated by squares, triangles and circles, 
respectively. It is implied that light-quark loops are inserted in all possible 
ways, for example as shown in the case of diagram (b). Also implied are 
mirror copies of the diagrams if appropriate. 
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transforms of these matrix elements contain ultraviolet renormalon poles at u = 1/2. For 
dimensional reasons, the residues are proportional to Ay times the matrix element p? G of 
the chromomagnetic operator, which is the only lower-dimensional operator contributing 
to the mass splitting. Therefore, the leading ultraviolet renormalon ambiguities of the 
p^ parameters are proportional to /igAm. Note that diagrams in which the external 
gluon is attached to a heavy-quark line, or to a two-gluon kinetic or chromomagnetic 
vertex, do not contribute. The diagram (cl) does not contribute to the chromomagnetic 
structure since the spin-orbit vertex only couples to the chromoelectric field. Moreover, 
with a suitable choice of the external momenta one can achieve that the diagrams (c2) 
and (c3) vanish. (Their sum vanishes with any choice of momenta.) A subtlety that 
needs to be taken into account is that there exists an ultraviolet renormalon at u — 1/2 
in the 1/m suppressed contribution to the quark wave-function renormalization with a 
kinetic-energy insertion, shown in diagrams (dl) and (d2); the corresponding ambiguity 
AZq = |Am/m in the normalization of p? G contributes to Ap\ G and cancels the Cp 
terms in the diagrams (al) and (a2). From a direct evaluation of the non- vanishing 
diagrams, using the method described above, we obtain for the ultraviolet renormalon 
ambiguities in the large-/5 limit Apf = —ki[C A j C p)p? G Am with coefficients k nG = 2/3, 
k A = 19/12, and kps = 1/2. The sum of the ultraviolet ambiguities in (0) is thus given 
by 

7C 

Apia + A PA - ApL = pi Am + 0(1/$) . (30) 

It precisely cancels the infrared renormalon ambiguity in the leading term. 

It is interesting that the requirement of a cancelation of renormalon ambiguities in 
physical quantities allows us to derive further information about the asymptotic be- 
haviour of the expansion coefficients c n in (g|) without any additional dynamical in- 
put |25|, [2~6fl . The point is that, beyond the large-/?o limit, the different terms in (Q) con- 



tain different powers of [a s (m)] 70//2/3 °, since they contain different powers of the Wilson 
coefficient C(m). These leading logarithms are exactly known to all orders of perturba- 
tion theory. They multiply the renormalon ambiguities of the various terms. In order to 
maintain the cancelation between infrared and ultraviolet renormalon ambiguities in the 
presence of the leading logarithms, the Borel transform for the Wilson coefficient in (f27|) 
has to be modified. In the vicinity of u = 1/2, the simple pole in 

S{u) = (2C F + ~C A ) I A- + . . . (31) 
must be replaced by a sum of branch points, 



u) i+/Wo 2 



?r k l r k i 19 CaK * i 1 CaK 



3 12 fI_ M )-W2/3o 2 (I_ U )W2A) 



(32) 

where 70 = 2Ca is the one-loop coefficient of the anomalous dimension of the chromo- 
magnetic operator given in (|4]), and Pi denotes the two- loop coefficient of the /^-function. 
The normalization constants K j = 1 + 0(1/ Po) are undetermined beyond the large-/?o 
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limit Up to corrections of order 1/n, the corresponding asymptotic behaviour of 

the expansion coefficients c n is 



c n+l 



(2(5 ) n n\n^^ 



2 IQ 
4C F Kt - -C A K 2 + -C A K 3 + C A K 4 n 

3 6 



7o/2/3o 



(33) 



For very large n, the last term gives the dominant behaviour (since 7o//3o > 0); however, 
for moderate values of n all contributions are of similar importance. 

In order to evaluate the Borel integral (|T5|), it is useful to rewrite it in the form ||T4" 



c(m) 



1 

A) 



<i In r ■ 



u>(r) 



/v w(r) 



lnr + In (m/Ay 
a s (v^e- 5 / 6 



O(l//? 2 ) 



47T 



(34) 



The function 

w(r) 
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2tH 



uq+ioo 



J du S(u) t u = C f w f (t) + C a w a (t) ; < u < 



2 ' 



(35) 



UQ— tOO 



which is the inverse Mellin transform of the Borel image S(u), describes the distribution 
of gluon virtualities in the one-loop diagrams in Fig. [I], which contribute to the calculation 
of c(m) . Using the methods developed in [[14]] , we obtain 



w f (t) 



W A (T 



- 4r + 2r 2 



1 + 4/r 
5r 

T 



1 + 4/1 



1 + 4/r - 1 



0(r 



(36) 



2^1 + 4/r 

These functions are shown in Fig. ^. Note that the integral in (plj) runs over the Landau 
pole in the running coupling constant, located at tl = (Ay/m) 2 . This is how infrared 
renormalons make their appearance. As in the case of the original Borel integral, we 
must specify an arbitrary prescription of how to deal with the Landau singularity. The 
renormalon ambiguity is given by the residue of the pole, i.e. A[c(m)] = w{tl)/ (3q. From 
an expansion of the distribution function in the region r < 1, we readily recover the 



result @S 



The representation fl3"4"[ ) makes explicit that perturbative calculations contain long- 
distance contributions from the region of low momenta in Feynman diagrams. At the 
same time, it provides for a convenient way to separate these long-distance contributions 
from the short- distance ones by introducing a hard separation scale A. We may then 
define the short- distance coefficient c s d(^, A) as 



Csd(m, A) 



A 2 



dy 
/i 2 



w(/j, 2 /m 2 



cue 



- 5 / 6 /i) 



47T 



0(1/ (3D . 



(37) 
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1.2 




Figure 3: Distribution functions u>f(t) (dashed) and Wa{t) (solid). The 
arrows show the factorization point for A = 1 GeV and m = 4.8 GeV (left), 
m = 1.44 GeV (right). 



As long as A is chosen large enough, this coefficient can be reliably calculated in per- 
turbation theory and is free of renormalon ambiguities. The long-distance contribution 
eliminated by this procedure must be combined with the non-perturbative power correc- 
tions in the heavy-quark expansion. 



4 Numerical results and conclusions 

Our goal was to study the large-order behaviour of the Wilson coefficient of the chro- 
momagnetic operator in the HQET Lagrangian. The renormalization-group invariant 
coefficient C(m) defined in (Bl) is known exactly at next-to-leading order ||. Here, we 
have derived an excellent next-to-next-to-leading order approximation by combining the 
exact two-loop matching condition obtained in || with an approximation for the three- 
loop anomalous dimension. This estimate has been obtained by summing fermion-loop 
contributions to all orders in perturbation theory, thus deriving the leading term in an 
expansion in powers of 1//3q. We have also derived an all-order result for the Wilson 
coefficient in the large-/?o limit, finding that the perturbation series for C(m) exhibits 
a divergent behaviour already in low orders, caused by a nearby infrared renormalon 
singularity. The resulting ambiguity is commensurate with terms of order 1/m 2 in the 
effective Lagrangian, whose corresponding ultraviolet renormalons we have identified. 

Let us now investigate the implications of our results by comparing the various ap- 
proximations for the Wilson coefficient. To this end, we evaluate the partial sums 



C( N )(m) = [a s (m)) 



70/2/3 



N 



n=l 



(38) 
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Table 1: Perturbative approximations for the Wilson coefficients C(N){ m b) an d C/m(m c ) 





C(N) 




v IY ) 




C( N ){m h )/C {N) {m c ) 


N 


large-/3 


exact 


large-/3 


exact 









0.580 


0.580 


0.692 


0.692 


0.838 


0.838 


1 


0.734 


0.659 


0.993 


0.847 


u. i oy 


fl 778 
U. 1 1 o 


2 


0.801 


0.697(1) 


1.207 


0.967(4) 


0.664 


0.721(3) 


3 


0.850 




1.463 




0.581 




4 


0.898 




1.875 




0.479 




5 


0.958 




2.714 




0.353 




Borel sum 


0.839(38) 




1.027(136) 




n SI 7+0-082 
U- 81 ' -0.065 




min. term 


0.898(48) 




1.207(214) 




n 744+0- 112 




fi > 1 GeV 


0.774 




0.906 




0.854 





for rif = 4 light quark flavours, using a s {rrib) = 0.22 and a s (m c ) = 0.36 in the MS scheme. 
The ratio of the two coefficient functions at these scales (for rif = 4) gives the perturba- 
tive correction to the ratio of mass splitting between the ground-state pseudoscalar and 
vector mesons in the bottom and charm systems ||: 

^ 0.89, (39) 

where A e g is a combination of the hadronic parameters pf introduced in (|7|), and for 
simplicity we neglect short-distance corrections in the 1/m terms. 

Our predictions for the Wilson coefficients C^){m b ) and C^){m c ) are given in Ta- 
ble |l| In the upper part of the Table, we show the partial sums obtained at different 
orders in perturbation theory, both in the large-/?o limit and exactly, so far as the exact 
results are known. In the large-/?o limit, the divergent behaviour of the perturbation 
series sets in already in low orders. At the scale nib, the minimal term in the series 
is reached around = 4, whereas it is reached around A^ = 2 at the scale m c . As a 
consequence, the values of the coefficient functions at the two scales become more and 
more different for larger N, and their ratio drifts away from the experimental value of 
the ratio of mass splittings. 

The reason for this behaviour lies in the nearby location of the first infrared renor- 
malon. This can be understood intuitively by considering the distribution functions 
shown in Fig. |3|, which are broad and extend far into the infrared region. Note that 



Mb* - M% _ U[m b ) 
Ml. - Ml = C(m.„) 



l n 



1 + A eff ( 



rrir rrih 
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the point /i = m corresponds to lnr = 5/3 in the MS scheme. Hence, there are es- 
sentially no contributions from scales above m (in the case of the Ca terms those are 
subtracted by MS renormalization), but on the other hand the distribution functions 
decrease very slowly (~ y/r) in the infrared region r — ► 0. The strength of this decrease 
is determined by the position of the nearest infrared renormalon singularity. In general, 
if the nearest infrared renormalon is located at u — k, then w(t) ~ r fc for r — > [ |i~4|j . 



As a consequence, the dominant contributions come from scales significantly below the 
heavy-quark mass.Q In the lower portion of Table [I], we show the results for the Wilson 
coefficients obtained by taking the principal value of the Borel integral. The quoted er- 
rors reflect the renormalon ambiguities, which are sizeable, in particular, at the scale m c . 
The Borel resummation gives values close to those obtained when the series is truncated 
at the minimal term (in the latter case the error is taken to be the size of the minimal 
term), as it is expected on general grounds. In the last row we show the short-distance 
contributions to the coefficient functions arising from virtualities above 1 GeV (in the 
V scheme), as defined in (|37|) . The portion of infrared contributions is much larger at 
the scale m c than it is at the scale mj. This is clearly seen in Fig. |3|, where the arrows 
indicate the location of the separation scale A = 1 GeV (corresponding to 435 MeV in 
the MS scheme). If the infrared contributions are cut away, the resulting values of the 
coefficients are again close to their values obtained from low-order calculations. 

Although the all-order results are instructive, we must not forget that they are ob- 
tained in a very questionable approximation scheme (the large-/?o limit). Indeed, the few 
exact results available indicate that the perturbation series may be much better behaved 
than indicated by the large-/?o limit. Still, we trust in the qualitative observations that 
the series start to diverge at some (low) order, and that the onset of the divergence is 
reached early the lower the scale m is. We will now explore the implications of these 
results for the phenomenology of the heavy-quark expansion. For a consistent inclusion 
of power-suppressed effects the perturbative coefficients of the leading terms must be 
known with sufficient accuracy. Different truncation or resummation schemes imposed 
on the short- distance coefficients imply different definitions of the hadronic parameters 
appearing at higher order in the heavy-quark expansion, such as the parameter A e g in 
(|39l). Only when the perturbative coefficients are truncated close to their asymptotic 
value, the hadronic parameters are m independent as they should. When applied to the 
particular case of the chromomagnetic operator, the conclusion is that there is not much 
to be gained by calculating C(m) beyond the first few orders of perturbation theory. 
Taking into account that the series for C(m c ) diverges earlier than that for C(rrib), we 
may argue that the optimal perturbative prediction for the ratio based on exact infor- 
mation is obtained by combining the NNLO result at the bottom scale with the NLO 



2 For the combined distribution function w(t) the average value of lnr is —1.310, corresponding to 
the scale /iblm ~ 0.226to in the MS scheme. This is precisely the scale obtained in the BLM scale- 
setting prescription [jlTj . In view of the low value of the average virtuality, the bad convergence of the 
perturbation series is not surprising. 
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result at the charm scale. This gives 

C(2){m b ) 
<5(i)( m c) 



0.80 . (40) 



In the large-/3 limit, the corresponding ratio equals 0.81 and is indeed very close to 
the principal value of the Borel integral. Using the result (|40f) , we conclude that the 
power corrections in (|39| ) give a contribution of about 11%, corresponding to a scale 
A c fj « 220 MeV. Thus, if the asymptotic behaviour of perturbation theory is carefully 
taken into account, it appears that power corrections in the heavy-quark expansion are 
well under control. 
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